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I 

^ ■ We discuss the multipolar expansion of the electromagnetic field with an emphasis on the radiated 

1—5 ; 

. field. We investigate if the employment of Jefimenko's equations brings a new insight into the 

calculation of the radiation field. We show that the affirmation is valid if one finds an interesting 

example in which inverting the order between spatial derivatives and integration is not allowed. 

^ ■ Further, we consider the generalization of the multipolar expansion of the power radiated by a 

Q ! confined system of charges and currents to a higher arbitrary order. 
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I. INTRODUCTION 



Despite the successful and long history of the electromagnetic field theory, there are 
several topics open to new theoretical and pedagogical contributions. One of them concerns 
the formalism of multipole expansion of the field, in general, and of the radiated one, in 
particular. Another issue is related to the importance of Jefimenko's equations in the study 
of such problems. Motivated by recent publications on the topic, we discuss some features 
of this type of problems. 

The multipole expansion of the electromagnetic field in Cartesian coordinates is exposed 
in electrodynamics textbooks, as the well-known Refs.- and^. Ordinarily, these expansions 
are calculated only in the first two or three orders, the higher-orders being considered too 
complicated. As Jackson writes in his textbook, the labor involved in manipulating terms 
in the expansion of the vector potential becomes increasingly prohibitive as the expansion is 
extended beyond the electric quadrupole terms (see Ref.-^, pp 415-416). For this reason and 
due to the applicability only in the long-wavelength range, another treatment, based on the 
spherical tensors and the solutions of Helmholtz equation is preferred. This alternative has 
also a larger domain of applications. Actually, starting from the results obtained employing 
this calculation technique, the reader can verify what effort is involved when returning to 
the multipole Cartesian moments which offer a higher physical transparency (see Ref.-). A 
relatively recent textbook^ and a paper-, the last related to the importance of Jefimenko's 
equations for expressing the electric and magnetic field when discussing the radiation theory, 
brought our attention on a very hard formalism employed for the calculation of even the 
first three or four terms of the expansion series. Though there are some prescriptions in the 
literature^,-,- for calculating higher-order terms of the multipole series based on a simple 
algebraic formalism of tensorial analysis, it seems that there is some reticence in using this 
last technique. For this reason, one of the aims of the present paper is to show how one 
can hide, as much as possible, the higher-order tensors behind some vectors, reducing the 
calculation technique to the formalism of an ordinary vectorial algebra or analysis. 

Another aim of the paper is to investigate if the use of Jefimenko's equations brings, 
as sometimes presented in the literature, a new insight in the calculation of the radiation 
field. We show that unless one finds an example where the spatial derivative and the 
integral operations can not be inverted, it is not always a necessary complication. As already 
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mentioned, we also use the opportunity to explain some advanced features of the multipolar 
expansions in the field radiation theory. 

We start in section [TTl by shortly presenting the notation convention we use and by giving 
a general formalism for handling multipolar expansions in Cartesian coordinates. In section 
llllt we derive the radiated electric and magnetic field without using the retarded potentials, 
while in section |IV] we present characteristics of the calculation for the radiation field when 
employing Jefimenko's equations. The advantages and disadvantages of different approaches 
are analyzed. In section |Vl we further discuss some features of the calculation for the radiated 
power, with an emphasis on the 4 — th order approximation in d/X < 1. Finally, in section 
\Yl\ we give the guidelines for the general tensorial calculus of the electric and magnetic 
moments. The last section is reserved for conclusions. 



II. GENERAL FORMALISM 



We write Maxwell equations with a notation independent of the unit system ("system 
free" Maxwell equations): 

„ „ jJo f ^ dE\ „ ^ IdB 



a \ dt J ^ a dt ^ 

V ■ B = 0, V ■ E = —p , (1) 

^0 

where Eq, /io, « are proportional factors depending on the system of units and are satisfying 
the equation 



a 
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(2) 

c is the vacuum light speed. Maxwell equations written in SI units are obtained from 
equations ([1]) for a = 1 and the SI values of Eq, /iq. For the Gauss system of units, 
a = c, Eq = l/47r, /io = 47r. With this notation, Jefimenko's equations can be written as^ 

/io /■„ [J] ,3 / 



V 



Eir, = --i- / vM Sx' - A / ^ d'.', (4) 

^ ^ Atteq J R Aira^ J R ^ ^ 

V V 
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where R = r — r', [/] = f{r',t — R/c) and the origin O of Cartesian coordinates is in the 
domain V. The support of charge and current distribution is supposed included in T>. If in 
equations and (jl]) the order of the derivative and the integral is inverted, one obtains 
the well-known relations between fields and potentials: 

1 (94 

B = VxA, E = -V<^--— , (5) 

a at 

with the retarded potentials 

A(r, t) = ^ [^4 $(r, t) = [ 4 dV . (6) 

^ ^ Ana J R ^ ^ Attsq J R ^ ^ 

V V 

In^, the authors derive the multipole expansion of the radiation field from equations 
([3]) and (jlj) claiming to give an original demonstration specific for Jefimenko's equations, 
without employing the retarded potentials. We should agree with this claim if at least 
some calculation of the authors is different from those employing the potential multipole 
expansions which are generally used in literature. In the following, we search for a difference 
between the calculation presented in Ref.-^ and the standard one making use of potentials. 
The goal of the exposition below is to inform on some results regarding multipolar expansion 
in Cartesian coordinates, too. 

Let us derive the multipolar expansion of the field B given by equation (JSj), and written 
explicitly as: 

/^O „ , /" rj [Jk] -,3 / 



B{r, t) = / d, ^ d'x', (7) 



where e, are the unit vectors of the Cartesian axes. Writing the integral from equation 



as 



5,I#dV=/5,(Ml^| d^x', (8) 



R J \ R 

we obtain the multipolar expansion of the magnetic field about O as function of r' using 
the Taylor series of the integrand: 

n>0 ■ ^ ^ 

Equation ([7]) can now be expressed as 

- n>0 ■ ^ ^ 



where we employed the notation [/]o = f{r', t — r/c). 

We assume one is allowed to invert orders of operations in equation ([9]) and to write: 

= ^ ^-^^^^'^ • • • J ^ / • • • < {■J'^]o ) dV. (10) 



Equation (fTOj) represents the curl of the multipolar expansion of the vector potential A. 
Thus one can perform firstly the multipolar expansion of this potential. It is the usual 
procedure. 

No matter what procedure is employed, a constant in the calculation is the presence of a 
vector a(r, t; n) defined by the Cartesian components: 

ak{r, t- C, n)=Q,... C„ ^ j . . . < Mr', t) d'x'^ . (11) 

Here, ( can be either an operator or a number. Generalizing to the dynamic case an algorithm 
used in- for the magnetostatic field, we introduce in equation flTTl) the consequence of the 
continuity equation, written for to = t — r / c: 

Jk{r',to) = V'{x'j,J{r',to)) +x'J{r',to). 

We obtain 

afc(r,to;C^) = Cn---0„^y" < . • . < V (x'JJ]o) dV + Cn • • • 0„^Pn....„ fc(to), (12) 

V 

where the Cartesian components of the n — th electric moment of the given charge distribu- 
tion: 

PH...^At) = j <...<p(r',t)dV (13) 

V 

are introduced. Let us define, for simplifying the notation, the vector 

P(r, t- C, n) =e,Cn... ^^^^i^. (14) 
Performing partial integration and taking into account that J vanishes on the surface dT>, 
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equation (fT2i) can be written and processed as follows: 

afc(r, to; C, n) - Vk{r, t; C, ^ + 1) = "Oi • • • ^^-^ / ['^]o ' ^'^^n • • • O d'^' 

= C. . . . C..i / X- . . . xL. (4 l^do - <. 141.) d'x' - na.(r. «„; C. n) 

= -neki„gCinCh---Cin-i^ J < •••a;-„_^(r' X [J]o)^dV-nafc(r,to;C,«)- (15) 

V 

Introducing the n — th order magnetic moment, as in Ref.'^, by its Cartesian components 

/ (r'xJ(r',t)),^dV, (16) 

V 



Mn...i„(t) 



{n + l)a; 
equation (|T5l) becomes 

afc(r,to;C,^) = -«^fci„qCi„Cii---0„-i-'^^^^^^^^^^^^ + ^-T^fc(^'^o;C,^+l)- (17) 

r n + 1 



Similarly to equation flMj) . we introduce the vector 



M(r, t; C, n) = e, Cn ■ ■ ■ , (18) 



writing finally equation f[T7j) as 

a(r, to; C, ^) = -«C X M{r, to; C, n) + ;^'^(^> ^o; C, « + 1) • (19) 

With this result, the magnetic field from equation (fTOj) can be expressed with the help of 
the vectors Al and 7-*: 



B{r, t) = -^^f^i^ijk afe(r, to; V, ra) 

n>0 

- £ E X (V X M(r, U, V, n)) + £^ E (^T^V x V(r, V, „ . 1), 

n>l n>0 ^ ^ 

or, by a change of the summation index in the second sum, 

B{r, t) = V X ^ V tlLl ("v X A4(r, to; V, n) + -P(r, to; V, n)] . (20) 

n>l ^ ' 



6 



From the last expression one has no problem in identifying the multipolar expansion of the 
vector potential A: 

Air, t) = ^Yl ^^n^ ^ ^' ^) + ^' ^)') • (21) 

The calculation for the electric field can be performed in a similar manner. One obtains: 

^ n>0 

-7^ E ^^n— ^ ^) + ^' (22) 

47ra; ^-^ n\ \ a I 

n>l ^ ' 

Comparing equations ([5]) with (1211) and (l22l) . we single out the multipole expansion of the 
potential 

$(r, t) = 5^ ^-^ V ■ P(r, to; V, n). (23) 



III. RADIATION FIELD 



For calculating the radiation field it is sufficient retaining only terms of order 1/r and 
for r oo. In most textbooks one retains only the terms of order 1/r, the goal being, 
usually, only the derivation of the radiated energy or of the linear momentum. Actually, 
when the goal is the complete definition of the radiation field, one must be able to derive all 
transferring properties, including the angular momentum loss. These are, in fact, minimal 
conditions for defining a physical system. In the last case, the terms of order 1/r^ are also 
necessary (see Ref.- -Problem 2 at the end of Section 72, and also Refs.-,— ). Although the 
aim of the present paper is different, we also give the formula for introducing terms of order 
1/r^ required for the evaluation of the angular momentum loss. The terms of the orders 1/r 
and 1/r^ are selected making use of formula^": 

'-JlX ■ ■ ■ '-'in I I „ '^n • • • '^in 



(^" • • • - '^{ii...i„_2'^in-ii„}) ^ g/-^^ • (2^) 

Again = t — r/c and z/j = Xi/r. By Ajj^ . j^} we understand the sum over all the 
permutations of the symbols iq that give distinct terms. The coefficients Dn are defined by 
the recurrence relations: 

= D^_, + n,Do = 0. (25) 



The formula from equation (12^ can be easily proven by recurrence. 

In fact, formula flMl) represents the sum of the terms corresponding to / = and 1 of a 
general formula 

a..-9<,.^ = i:4"lM;^^/W- (26) 

1=0 

In this equation, cl^'^l^{u) are symmetric coefficients expressed as linear combinations of 
products of components Ui and Kronecker symbols where i, ig, ip = ii . . .in- 

Considering equations ( |T4l) and ( |T8l) . one can see that P(r, to; C) and Al(r, to; C; ''^) 
are solutions of the homogeneous wave equation for r 7^ 0. Consequently, one can apply the 
formula for these quantities. Let us consider the multiple derivative as, for example, 
. . . , ^ , (-1)"+"^-^ M,,...,„_,fc(to) 

/ 1 \n+m— 1 Qj.n+m— 1 

+ ^( Vr^) = ■ ■ ■ -.^Q^r-^ir^ ^0; v. n) + (9(l/r^) (27) 

and similarly for T'{r,to;'V ,n). Using equation ( 1271) and the equivalent relation for 
{r , to] 'V , n) in equation (1201) . we obtain the first approximation of the multipolar ex- 
pansion of radiated magnetic field, which is sufficient for calculating the radiated energy 
and the linear momentum: 

1 a'^+i 



Bra,{r, = £ E 9^ X (^ X M{r, to; ^, n)) - x P(r, to; ^, n)) (28) 

n>l 

n>l 

For the field E^ad, we obtain 



1 ■ 1 gn+i 
^^'^'''^^ = ^ign+i • ^(^' ^0; i^, n) - P(r, to; n) 

^ n>l 

(y. 

+-U X M{r, to; z^, n)) . (29) 



We made use of equation ([2]). From equations (I2T1) and (1231) . using equation (l24l) . we get 
easily the expansions of the radiation field potentials: 

Aradir, t) = ^Y.^d^ (-^(^' ^, n) X ^ + ^P(r, to; ^, n)) (30) 

n>l 

and 

1 19" 

*-('-.*) = si;E;^s^'''''('-.*«;''.")- (31) 



Given the above expressions, one can verify that the relations between fields and potentials 
are: 

1 • 1 • c 

Brad = -{Arad X v) , Erad = - (Arad X v) X V = -Brad X V. (32) 

c a a 

These parts (proportional to 1/r) from the radiated electric and magnetic fields are purely 
transverse fields, satisfying the properties (see alsoi^): 

^■Erad = ^i ^■Brad = 0, £o\Erad\'^ = |-BradP- (33) 

IV. IS REALLY THE RADIATION FIELD CALCULATION FROM JEFI- 
MENKO'S EQUATIONS A NEW INSIGHT IN THE RADIATION THEORY ? 

In Ref.-, calculating the radiation field (in the first approximation), the retarded poten- 
tials are approximated by the formulae (66.1) and (66.2) of this reference, written here with 
the "system free" notation: 

$(r, t) ^ / p(r', to + -iy- r') dV, A(r, t) ^ / J(r', to + -u- r') d'x'. (34) 

AttEoT J c Airar J c 

V V 

As suggested in Ref.- (page 184, first footnote), introducing these expressions of the radiated 
potentials and retaining only the terms of order 1/r, one obtains the relations from equation 
( !32|) ( i.e. equation (66.3) from Ref.-). In this calculation, the derivative operators must be 
introduced in the integral and so, the proof is indeed realized directly for the fields E and 
B. For B, using the relations 

J{r'ta + ^) , l_, vt' 



and 



r r c 



V- = 0(1), V{to + -u-r') = --iy + 0{-), (35) 

r c c r 



we can write 



/io J{r',to + ^). , 



B{r,t) = Vx ^ ' " ^ ^ dV 

Ana ' 



V 



^° ^ux [ jir',to + -iyr')d'x' + Oih. (36) 
J c 

V 



Anac r 



From this last equation we can extract the term of order 1/r which represents the radiated 
field (precisely, the first approximation of this field) : 



(37) 



Bradir, t) = -— ^ i^x [ j{r\ U + r') A^'x'. 
Anacr J c 

V 



This is, in fact, equation (28) from Ref.- and we can consider, by examining Ref.", that this 
result was obtained long time ago from the expression ([3]) of the field B (see Ref.— and 
Ref.— ). Look also in Ref.- to see the usage of Jefimenko's equation for calculating -Brad- 
Introducing also the approximate expression E starting from equation (jll), we obtain 

E{r,t) = --^ / Uvp{r'M + -u-r') + J{r\t, + -u-r')] dV + 0(1) 
47ra^ r J \ c c J 

V 

V 

Using the continuity equation written in the point r' at the retarded time t — R/c, 

[p] = -[V- J(r', t')]t'=*-ii/c = - V • [J] + [ j] ■ V'(t - R/c), (39) 

it results 

p(r', to + -u- r') = - V ■ J(r', to + -i^- r') + ■ J{r\ to + -jy- r') + 0{-). (40) 

c c c c r 

The first term from the right-hand side of the last equation gives no contribution to the 
integral from equation (l38l) and, after a simple algebraic calculation, we obtain the expression 
of Erad from equation fl32|) (see also equations (66.3) from Ref.- and equation(30) from Ref.-). 

We point out that equation (!39|) or, generally, the relation between the space derivative 
of a retarded quantity and the retarded value of the space derivative of the same quantity, 
should be well-known for each student from a class of electrodynamics since when writing 
the retarded potentials as solutions of the wave equation, it is necessary to verify the Lorenz 
condition. The verification can be realized in a direct calculation, a good exercise for a 
student. Only in this way one can be convinced that the retarded solutions are indeed 
electromagnetic potentials. For this goal, equation fl39|) is indispensable since the Lorenz 
condition appears as a consequence of the continuity equation. 
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Based on the above results, we can quickly obtain the multipole expansion of the radiation 
field. Considering the adequate Taylor series for the integrand in equation (|36|) . we write 



X E • • ^ /< • • • <J(r\ to) dV 

n>0 



= -^1^ ^dfTT ^ X «(^' (41) 

n>0 

The vector a was defined in equation (fTT]) . After introducing its expression as a function 
of the vectors associated to the electric and magnetic moment, equation (fT9|) . the radiated 
magnetic field becomes 

n>l 

n>0 ^ ' 

The final step is changing — > n — 1 in the second sum of equation We obtain again 
equation fl28l) . 

We can say that the result fH2|) . which is the multipole expansion of the radiation field, 
was determined via Jefimenko's equations. 

Alternatively, after reaching the expression given by equation ([9]), we can perform the 
multipolar expansion of the radiation field in such a manner that we can also say that it 
is obtained via Jefimenko's equations. We extract from the integrand in equation ([H]) the 
terms of the order 1/r (and for a complete definition of this field): 

Bradir,t) = -^e,e,,ki^,v,,...Vi^- / ^ a;:^ . . . x^^^^fj^jo dV 

^ n>0 

^0 1 



n>0 



arriving, as expected, to equation ( 14T1) . 

We remind the reader that for obtaining the results of equations (1201) and ([22]) we admitted 
the commutation of the derivative with respect to the spatial coordinates, with the series 
expansion and the integral on the domain V. In the present section, for the case of the 
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radiation field, we avoid the inversion of tlie derivative witli tlie integral operation. Regarding 
the commutation with the Taylor expansion, we consider that such an operation cannot be 
avoided as long as we want to emphasize the multipolar moments. 

As a short extension in the argumentation of the utility of Jefimenko's equations, we 
show how one can avoid the commutation between derivative and integral for an arbitrary 
point in the exterior of the domain V. Indeed, using equation (l26ll in equation ([9]), we can 
write the magnetic field with the help of time derivatives 

n>0 ■ /=0 

Here, 

V 

Similar to the derivation of equation (JTTj) from equation (ITT]) , we get: 



n 



...in k' 



2 J ll...tn n---i7i 



Inserting this equation in equation fj44l) . we obtain 



n>l ■ «=0 

_^ /^O ^ (-1)" _ ^ 1 



in + i)!^^^^^'^Z^ r^+i^^^-^i--" 5^n+2-;^^i--nfc ^^^^ 

n>0 ^ ' ^ Z=o 

As anticipated, it represents the multipolar expansion of the magnetic field for any point in 
the exterior of the domain T) and it is obtained without inverting the spatial derivative with 
the integral operation. For the approximation of the radiated field (in 1/r), one can easily 
verify that writing equation fj45l) for / = and taking 

we get equation (!28l) . 

In conclusion, employing Jefimenko's equation in the radiation theory could bring a new 
insight only if the inversion of the spatial derivative and the integral operation is not allowed. 
We admit we were unable to find an interesting example where an inversion is not permitted, 
at least for generalized distributions. However, it might be possible to find such examples, 
and, in this case, the indispensable character of Jefimenko's equations would be obvious. 
Otherwise, for the regular cases, it appears as an unnecessary complication. 
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V. SOME FEATURES OF THE RADIATED POWER CALCULATION 



The Poynting vector is 

S=—{ExB) = —\Brad\^ u + Oil/r') = eoc\Erad\^ u + Oil/r^). (46) 

The total radiated power may be written as the hmit of a surface integral on a sphere 
centered in O, of radius r, for r — > oo. Let us express the energy current in the radiation 
approximation corresponding to the sphere of radius r at the time t: 

iV(5,E,;t)= (( r^u- Sir, t) dn{u) = — ((r^\Bradir,t)\^dn{u) + 0{-). (47) 
J 1^0 J r 

The quantity A^(<S', S,.; t) dt represents, therefore, the energy which crosses the sphere in 
the time interval (t, t + dt) and is determined by the values of the multipolar moments of 
the source in the interval (t — r/c, t + dt — r/c). For large but finite r, the integrand from 
equation ( l47l) can be enployed for drawing conclusions on the electric charge distribution at 
the retarded time from observations on the angular distribution of radiation. Since dt = dto, 
we can say that A^(>S', S,,; t) dt is the part of the energy emitted by source in the given time 
interval which contributes to the energy intensity corresponding to at the time t. As 
one can see from equation (l28!l . if in equation (H7|) we put t instead to, then the quantity 
lim A^(>S', S.^; t) dt represents that part of the energy emitted by the source which contributes 

1 — ►oo 

to the radiated energy or, shortly, radiated by the source. The situation changes when the 
support of the source depends on time and, in particular, for the radiation of a moving point- 
like source (seel, & 73). in conclusion, the energy emitted by the source is characterized by 
the intensity 

1= lim N{S,K;t) = — / \Brad{r,t)\^dn{u) = {r''\Brad\^), (48) 

where 

Bradix, t) is given by equation ( 128|) substituting the retarded time to by t. Since in equation 
fHHj) the factor is simplified by the factors 1/r included in the definition of M. and V we 
should rather introduce the quantities ^ and tt: 

IJi{t]u,n) = r M.{r,t]u,n), n{t;v>,n) = rV{r,t;v>,n). (49) 

13 



The expression for the radiation intensity T is 

^W = 4^ 2. c-^- ar+i \ • ^ ~ ^, - - . X 7r(t; n 

n,m>l 

^m+1 , C \ 

( (i^ ■ i^, "^)) - i^, m) - - X 7r(t; z^, m)^ . (50) 

Employing the notation = 5"/(. . . ,t)/dt^, and specifying only the argument n in /x and 
TT when there is not of confusion, we can write 

, sr-^ 1 

n,,m>l 

(m)) 

+ ^,n+l(^) ■ A^,m+l("^) + -^J',n+lin) " (l^ X 7r,„+i(m)) 

c 

+ - (i^ X 7r,„+i(r2)) ■ /x„+i(m) + — {7r,„+i(n) ■ 7r,m+i(m) 

-(i/ ■ 7r„+i(n)) (i^ ■ 7r,„+i(m))}) . (51) 

One can calculate the averaged quantities from the last equation using formula^: 

/ \ / °' n = 2k + l, 

{^h---^ij=< ^ _ _ (52) 

I (2fc + l)!! '^{«l«2 • • • ^in-lin}} ~ A; = 0, 1, . . . 

In equation (15T1) all terms containing an odd number of factors v vanish and we can retain 
only terms with an even number of these factors. We point out that ^{n) or 7r(n) contain 
n — 1 factors z/ (see equations ([HD, (ITH]) and fll^ ). 
We have to calculate expressions as, for example. 



•Jm-l J- 



This is cumbersome even for the first approximations. 

Before discussing and applying the above results, we stress an essential issue for the 
existence of a precise approximation criteria when a finite number of terms in equation (182|) 
is retained. Let us consider the source from V being a system of point-like electric charges 
qi . . .qn- Therefore, 

N N 

p{r,t) =Y^q,6{r~ r»(t)), J{r,t) = g,f»(t)5(r - r»(t)), (53) 



i=l 4=1 
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where r^^\t) represents the position vector of the particle i. In this case, 

N 

P^....^.(t)=J2^^xf;(t)...x'{;(t) (54) 
i=l 

and 

N 

Mn...„W = E^^^n W • ■■^t.it)i^^^ X (55) 

1=1 

Let us suppose the particles oscillating with a pulsation u = 27ic/X, i.e. 

r»(t) = r^^^e'"*, f«(t) = iu r^^^^e'^* = ^ r^^e''^'. (56) 

A 

(For the general case of N wave lengths Aj, i = 1 . . . A^, in the following, we will understand 
by A the shortest of them.) Denoting by d the linear dimension of the domain T>, we have 
|*"o^| ^ d. Considering the radiation in the case of a long wave-length, A > c/, we introduce 
the parameter 

C = ^<1. (57) 
For the amplitudes of the source, we identify the orders of magnitude: 

^^C, |r«|~^<C. (58) 

Obviously, the same relations can be written for any x[''\ xf'\ For the time derivatives of 
the tensors P'-"'' and M we can conclude that 

With the approximation criteria considered above, we start by calculating the terms from 
equation (l5Ti) up to the A—th order in the parameter (. For the electric and magnetic dipolar 
moments we use the usual notation p and m, respectively. We select from equation flSTl) all 
nonvanishing terms for {n,m) = (1, 1): 

111 = ^5 {-{■Ui-Uj)mimj + fh? + — [p^- {i^ii^j)pipj) ) . (60) 

Inserting the result for the even combinations {viVj)., we obtain the well-known expression 
for the radiation in the dipolar approximation: 

a" 



15 



Now let us consider in equation (15T1) the terms with (n, m) = (1, 2) and (2, 1), discarding 
the vanishing ones : 

Xi2 +X21 = ■ (i^x 77 (2)) + {ux 7r(l))- A* (2)) . (62) 

The first term, written exphcitly with omission of dots representing the time derivatives of 
no significance for the tensorial relations, is 

■ (1/ X 7r(2))) = (mi(i/ x 7r(2)).) = eijkmi{vjVq)Pqk = ^eijkSjgniiPgk = , 

because of the symmetry of the electric quadrupole moment. For the second term, 

((i/ X 7r(l)) ■ Ai(2)) = {Eijki^jPki^gMqi) = ^pkSkijMji , 

and finally, 

J12 +X21 = ~~[^^:^^^~^Pk^kij Mij . (63) 

Writing equation flSTl) for [11,171) = (2,2), we discard the terms of order 6 in ( as, for 
example, /x (2)- ^ (2). In this approximation, 

^ PikPjk -{ViViPjUq) PuPq^j + 0{0 



1 /I 1 



1 2 

PijPij ~T7 PiiPjj ^TT P*iP*i ) "I" 



IGvreocS VS ' ' 15 15 



(PuPu -^(Pn)') +^(C') • (64) 



SOvreoC^ 

Terms of order 6 in C exist also for (m, n) = (1, 3) or (3, 1). Retaining only terms of order 
4, we get: 

2:13 + X3I = ^^(1) • 77(3) - {U . ^(1)) {U . 77(3)) ^ + O(C^) 

1 



{l^ll^qPi)Plgi - {l^iPjl^ll^q)PiPlqi ) + 0{C ) 
1 

?PiP«gi - Y^S{^j6lq}PiPlgj + C(C ) 
1 



P^ Pqqi + 0{C) . (65) 



QOvreoc 
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Adding up equations (16T]) . (163|) . (16^ and (1651) . we obtain the A — th order approximation of 
the total radiated power: 

1 /2 „ a ... 2 ■'■ 1 / 1 -2 \\ ^ ^ 

% = + - +iH?^^P- + ^ ^P^.P^. -3 P..J J (66) 

The output can be partially compared with a well-known result from literature (see Refs.- 
and^), but for this we have to introduce the irreducible electric and magnetic momenta 
defined as symmetric trace-free ( "STF" ) Cartesian tensors. Let us consider a n — th order 
tensor T*-"-* and the corresponding projections iS(T^"'') and T(T*-"'^) on the subspaces of 
symmetric and STF tensors. For the case of the electric moment P*-"-*, this is a symmetric 
tensor and one has just to establish their STF projection. Let us consider the simplest case 
of the quadrupolar electric moment P'-^-*. Writing the components Pjj as 



there is a unique value of the parameter A such that U^"^^ = T{P^'^^). For A = Pqq/3 , 

1 f 1 

Llij = Pij - -PqqSij = / {xiXj - -r'^Sij) pd^x . (67) 



V 



In equation (l66l) the octupolar electric moment P^^-* is present. The STF projection can be 
calculated searching the first order tensor A*-^-* such that the STF projection fl^^) = T(P^^'') 
is given by the components 

riijfc = Pijk — S{ij/\k} ■ (68) 
From the condition of vanishing traces of the tensor fl*^^-', one easily obtains: 



= -P,,. = - I r'xipd'x . (69) 



V 



Concerning the magnetic quadrupolar moment M*^^-', we have a simple procedure for STF 
projection. Let us write the identity 



Mi, = ^(M,, + M,,) + ^(M,,-M,0 



where the first bracket represents the symmetric part of this tensor, and the second one, 
the antisymmetric one. The symmetric part is, for this case (n = 2), a STF tensor T^^) = 
r(M(^)). Therefore, 

My = Vij + ^Eijk^k , (70) 
17 



where 
Since 



^fc.,M,,. = A J (rx(rx J)),d3x = A J {(r ■ J) r - j) ^ d'x . (71) 



V V 



we obtain 



P^jP*. pL=n,,n., . (72) 



Combining the magnetic quadrupolar and the electric octupolar terms from equation (l66l) . 
with their expressions from equations (jSZD, fpS]) - and (1701) . we get 

a ••• 2 •'■ 4 / a 1- \ 4 

a consequence of the traceless character of r\^^\ U^^^ and r^'^\ Here, the vector 

t = ^N-iA=l J {{r.J)r-2r'j) d'x (73) 

is introduced. The last expression is obtained from equations (17T1) . (l69l) and applying the 
continuity equation together with an operation of partial integration. This is the so-called 
electric toroidal dipole moment — ,— . We can write the radiation intensity X(4) in terms of 
STF projections of electromagnetic momenta: 

1 /2..2 2^2 .. 2 1 4 .. -A 

Now, we can compare this last result with the one given by equation (71,5) from Ref.- 
where the term corresponding to the quadrupolar electric moment is written in terms of 
the STF tensor D*^^) defined by the components Dij = Sflij. We see that the expression 
given by our equation ([71]) differs from the equation in Ref.i by the term — (4/3c^)p- t, i.e. 
the contribution of the toroidal electric moment. In Ref.-, one calculates the contribution 
to the radiation of electric and magnetic dipolar and electric quadrupolar momenta. For 
this calculation, one uses the STF tensor D^^^ based on the invariance of the field to the 
substitution Pjj — Djj or, equivalently, the gauge transformation of potentials for such a 
substitution. We point out that this invariance is a singular case and, for higher orders 
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of the multipolar expansion it is not true. The toroidal electric moment does not appear 
because the orders of magnitude of the different terms in the multipolar expansion are not 
consequently evaluated. The same omission is done in Ref.-. In Refs.— ,— ,-, the toroidal 
dipole is obtained as a first term from a class of toroidal multipoles. The class is pointed out 
in a complete multipole analysis, appearing as additional terms beside the STF projections 
of the primitive moments P and M. In Ref.—, our formula (17^ is obtained as a correction to 
the result of Ref.- in order to assure a quantity origin independent, i.e. translation invariant. 
The toroidal contribution from equation fl7^ is obtained in Ref.—, too, by a calculation 
close to the procedure from the present paper. 

Let us consider the contribution of the electric quadrupolar moment which, as seen from 
equation (1251) . is proportional to the third time derivative of the vector 

V = U X T = ^eiEijkl^jUgPgk- (75) 

The substitution 

p(2) ^ ^(2) ^ ^(p(2)) (76) 

in equation (175!) gives 

1 1 
V - —eieijki^jUg6qkPii = V - —Puiv x u) = V. 

Sr Sr 

One can see that B{r,t) and E{r,t) given by equations ( |20l) and ( !22l) are also invariant 
if one performs the substitution fITBl) . It is an exercise for the reader to prove that the 
substitution (17^ in the expansions of the potentials A and $ given by equations (j2ID and 
( !23|) has as effect a gauge transformation of these potentials^ 

Unfortunately, this type of invariance is not valid for magnetic moments and, for n > 3, 
for electric ones. In these cases, the physical results are not invariant with respect to the 
transformations 

pW^r(P")), M(")^T(M")). (77) 

Let us suppose however that, at least for a given pair of numbers (M, N), all the magnetic 
and electric moments for m < M, n < N in the expression (15T1) of the radiation intensity can 

~ (m) ~in) 

be substituted by STF tensors M , P • These new tensors do not necessarily coincide 
with the corresponding projections 'T(M*^™''), 'T(P^"'') as it will be seen in the following. In 
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this situation, the calculation of the terms from equation (15T|) obtained by the corresponding 
substitutions /x(m) Ji{m), Tz{n) tt (n) is more easily performed. Indeed, if A^*^-*, B^™'^ 
are STF tensors, then^S 



(7^ 



(2A; + 1)!! ~ (2A; + 1) 



Here, since A*^"^ and B^"*-* are symmetric tensors, it is of no importance what is the order 
of factors in the contraction A^^^.B^™^ ( noted as A^"^||B^™\ too). Obviously, the result of 
a such contraction is a STF tensor of \n — m\ — th order. All the terms in equation (IHTj) 
can be calculated using equation ( ITHj) except for those containing vectorial products which 
vanish as consequence of the equation fl52l) and of the symmetry properties of the tensors 
p(")^ g^ppiy equation ([78]) to the term : 



(2n - 1 
and 



( [u . nit; u, n)){u. J^it; u, m)) > = (m^"^ . M^^^^ (80) 

The expressions for alike terms are similar. The terms including vectorial products as, for 
example: 

{Ji{t]u,n) ■ [u X 7f(t;j/,m))) = {]li{t]h',n)eijki^jiTkit] 1^,171)) 

= £ijk {^h ■ ■ ■ ^i„-l ^jl ■ ■ ■ ^j„-l) Mil...i„-liPjl...j,n~l fc = (81) 

vanish, since each term from the average (i/j^ . . . i^j I'ji ■ ■ ■ T^jn^^) given by equation (!52|) 
includes either g = 1 . . . — 1 or g = 1 . . . m — 1 or directly an odd number of v. 
Using these results, we obtain the final expression of the (total) radiated power: 

where we consider the limit for M, N ^ 00. We point out that considering the infinite sum 

~{n) ~(n) 

in equation (1821) . makes the tensors M , P infinite series. This is not so catastrophic since 
for practically applications only a finite number of terms from equation (l82l) is necessary, 
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as it will be seen in the following. Obviously, one can apply equation (159!) for the STF- 
projections of the multipole momenta. 

For an answer to the question how one can use equation (1821) for calculating the radiation 
intensity, let us return to the problem of calculating this quantity up to the A — th order in 
the parameter (. Since, together with P*-^\ the momenta M*-^^ and P*-^-* give contributions 
in this order, we have to consider in equation (182!) the STF projections V^^^ and r\^^\ too. 
This time the field is not invariant for the substitutions 

M(2) ^ r(^), p(^) ^ n(^). (83) 

Let us calculate the effect of the first substitution from equation (155]) in equation fl25]) for 
expressing the radiated magnetic field. The terms from the expansion series of B affected 
by the substitution are 

1 1 1 

^ ' ^^'^ ^ - ■^('^' 2)} = {e^l^il^jl^l Mij - e,Uj M,-,} 

M{2)^r(2) 1 9^ 1 f 1 M r/, 1 

— ' 2(^dPr I ^'^'^^^^ ~ -^(^it^i^jt^ii^ijkf^k - eiUjMji + -BiUjSjikNk 

= ^ ■ ^^^'^ " ^^"^^ ^' 2) } - ^ ) ' 

where A?" = CjNj = From equation (I28p we can see that the alteration of B by the 

substitution M*-^^ f^^^ is compensated by the following transformation of the dipolar 
electric moment: 

Let us consider the substitution P^^^ fl^^^ in equation (1281) . The term affected is: 

P(3)^n(3) 119^. X ^ \ 

I , 11 



[u X P(r, to; 1^, 3)) + -—- 1/ X A 



where A = CjAj. The alteration of B by the substitution P^^-* fl*^^^ is compensated by the 
following transformation of the dipolar electric moment: 

p - P + ^ A. (85) 
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The total alteration of B by the two substitutions (183|) is compensated by the transformation 

4 6 / ^ c2 



P^P = p--J-N--\] =p--,t, (86) 



with t defined by equation fl73|) . 

~(i) _ ~(i) 

Now, in equation (1821) we have to consider P = p given by equation (1861) . M = 

~(2) ~ (2) ~(3) 

m, P = n*^^\ M = f^^^ P = n*^^^ and, for n > 4, the primitive momenta which 
contribute only with terms of orders greater than 4 in C,. We can write: 

f 2 / . 2 c2 1 ... 2\ 1 -(2) ...(2)) 

The quadrupolar magnetic and octupolar electric terms are not written because these are 

...2 

of order 6 in C,. Since t ~ C^, from the last expression of I we are left with equation (17^ . 

This calculation offers a suitable example for understanding the general scheme of ten- 
sor reduction for electric and magnetic moments, as well as for the usage of formula fl82|) 
when describing the radiated power. Replacing an electric moment of order n by its STF- 
projection induces in the tensors of inferior order, compensating terms of order C^. The 
general term in the expansion ( |82i) is of order (^^". If we are interested in the approximation 
of 1 up to terms of order C,"^^, then we have to replace the tensors related to the moments 
p(2fc-i) ^^^^ |y^(2fc-2) STF-projections and to take into account all compensating terms 
for the lower-order tensors. 



VI. SOME BASIC FORMULAE FOR GENERALIZING THE PROCEDURE OF 
REDUCING THE ELECTRIC AND MAGNETIC MOMENTUM TENSORS 

In this section we list the required formulae for the generalized tensor reducing procedure. 
The equations are given without the related proofs since they can be found in the literature. 

Let S'-"'' be a symmetric tensor of rank n. Its STF- projection results from the following 
formula: 

r(SW)^^ = S,,..,„ - <5|,,,A(S(")),3...,„}. (87) 

The operator A is defined on the account of a formula for the STF-projection of a symmetric 
tensor given in^ (the book— is cited as the origin of this formula): 

AfS(")) '"f^^' (-l)-[2n-l-2(m + l)]!! 

m=0 ^ '" 
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The proof can be found in^^. [a] stands for the integer part of a and the notation S^"'^^ 
indicates p pairs of contracted indices. In the following, for simplifying the notation, all ar- 
guments of the operator A should be considered as symmetric tensor i.e. A(T*-"'') = A(iS(T*-"^) 
for any tensor T"^ . The same applies to the operator T: r(T(")) = T(5(T(")), S being the 
symmetrization operator. 

In the symmetrization process we have to calculate the symmetric projections of some 
tensors L^"^ of the magnetic moment type: they are symmetric in the first n — 1 indices and 
the contraction of i„ with any index ig, q = 1 . . .n — 1 gives a null result. For the symmetric 
projection of such a tensor, we introduce the formula: 

^ n—l 

('5(L("^)).,,...„ = L.,...„ - - E^--X'1.„-..(L^"^ (89) 

A=l 

where Af\'^^ is the component with the i\ index suppressed. The operator Af defines 
a correspondence between L*-"^ and a tensor of rank (n — l) of the same type from the 
symmetry point of view. Particularly, 

^2k+l /,^(n)N _ {-^fn r 2^fc ^n-2k-l ^ X J) d^X, = 0, 1, 2, . . . (90) 

{n + l)aj 



V 

where a" x 6 is the tensor defined by the components (a" x b) . . = Oj, . . . (a x b)i^. 
Let us consider the process of STF reducing the multipole tensors starting with the rank 
+ 1 in the electric case, and in the magnetic one. In Ref.~ one finds a general formula 

~(n) ~(n) 

for the resulting tensors P , for n = 1, . . . , N and M , for n = 1, . . . , — 1. It includes 



all compensating terms obtained in this process 



[(N-n)/2\ 



\k Q2fc-1 



(n) 



P _ T(P(")) + ^ ^ 2 Q-f2k-l """fc ' 

k=l 

/ fc— 1 

-W _ (_^\k^2kn'i /iWAfc/^D^^+^'^^^i _^ 

1=0 



t(") = {-Ifc^^ri 4")a'=(P^"^''^) + ,A'Ar2'=-2i-i^,^(n+2/c-i)^ ^ ^9^) 



and 



M^"^ = r(M(")) + r |i5^c'i;^A'Ar2'=-2'(M("+2fc)) i . (92) 



k=l 1=0 
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The coefficients are given in a compact algebraic form in Ref.— : 



(n) 



1 




n{n + 2l)\ 



k,l 



(-1)^-' n(n + 2/)! 



(93) 



2«c2fc {n + 2k){n + 2k)\' 



The reader is encouraged to apply these formulae for the cases iV = 4, 5, 6 and to find 
convenient intermediary calculation. We point out that in equation ( 19T1) the normalisation 
of the quantities T*-"^ is chosen such that the quantities T^"'' coincide with the electric toroidal 
moments given in literature at least for up to 3. 

VII. CONCLUSION AND DISCUSSION 

After an introduction on the formalism for multipolar expansions of the electric and 
magnetic field, we dedicated sections IIIII and [IV] to the main purpose of this article: the 
analysis of different methods for calculating the radiated field with an emphasis on the 
novelty the use of Jefimenko's equations can bring. In the last two sections we presented 
aspects of the radiated power calculations from the tensorial point of view. 

From the analysis of sections IIIII and IIVI we can draw some conclusions on the utility 
of Jefimenko's equations when considering the multipolar expansion problem. As one can 
notice, when deriving equation (ITO!) from (Q, no matter if we work with potentials or di- 
rectly with fields, the inversion of the spatial derivative with the integral on the domain 
T> is mandatory. In Section IIVI it was proven that the multipole expansion of the fields E 
and B can be obtained generally and directly from Jefimenko's equations. All expanding 
operations are performed on the integrand, but, still under the assumption that the integra- 
tion operation is distributive with respect to the Taylor series. It remains only to argue the 
necessity of the corresponding additional calculation effort for applying this procedure. 

As an additional remark, we emphasize that in the same sections we tried to remind the 
reader that if one wants to completely describe the radiative systems, one has to include 
besides terms of order 1/r, the contributions from the field expansions. 

Ref.- is part of a paper series trying to emphasize the theoretical and practical importance 
of Jefimenko's equations. These equations are considered as extraordinarily powerful and 
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illuminating as the authors of Ref.— write. We have nothing against the open enthusiasm 
in these papers. We neither dispute the beauty of the result regarding the calculation of 
the retarded fields E and B directly from Maxwell's equations, without having to introduce 
and handle the potentials. Maybe a series of applications based on these equations are more 
efficacious and physically more transparent. Although, we are circumspect concerning the 
axiomatic treatment of the electromagnetic theory starting from these equations (opposite 
to opinions from e.g. Refs.— ,— ), but, this subject will be discussed elsewhere. 

In the second part of the paper we gave a detailed calculation for the radiated power 
up to the 4 — th order in d/X < 1. The method can be similarly applied to higher-orders. 
We related our procedure to existing prescriptions in the literature and we underlined what 
omissions might show when particular terms of the multipolar expansions are neglected. We 
hope we were able to convince the reader how powerful and not so complicated a consistent 
vectorial/tensorial computation is. 
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